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We revisit a weak gravitational lensing problem by constructing a setup which describes the actual 
system as accurately as possible and solving the null geodesic equations. Details are given for the 
case of a universe driven only by a cosmological constant, A, which confirm the conventional results: 
The conventional lensing analysis is correct as it is, without any need for correction of 0{h.). We 
also treat the cases of the lensing in generic Friedmann-Robertson- Walker backgrounds. 

PACS numbers: 95.30.Sf, 98.80.Es, 98.62.Sb 
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I. INTRODUCTION 

In the conventional analysis on the weak gravitational 
lensing, the bending of light from a source by a mas- 
sive deflector, or a lens, is calculated in the comoving 
frame with the Schwarzschild metric, whereas the effect 
of the propagation of light in an expanding universe is 
taken care of by using angular diameter distances as the 
distance measure. The underlying idea for this separa- 
tive treatment is that on most of its trajectory a photon 
does not know the existence of the lens, and it feels the 
gravitational pull only near a small region of the closest 
approach to the lens where the spacetime can be locally 
approximated by Schwarzschild. Although this sounds 
like a reasonable approximation, it is somewhat qualita- 
tive, and hard to produce quantitative errors to justify 
itself. 

Then a better, or at least more precise, description 
of the situation can be given by 1. putting our lens in 
a Friedmann-Robertson- Walker (FRW) spacetime to get 
the "Schwarzschild-FRW" spacetime and 2. solving for 
the null geodesies connecting a source to an observer in 
this background. 

Step 1 was done about 80 years ago by McVittie[l|], 
who got the exact solutions for the Einstein equation 
sourced by a localized (spherical) mass and cosmological 
medium. The McVittie solution is given by 
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where X is the usual comoving coordinate, 

m 

^~ 4a(i)|X-Xor 

with Xq the location of a mass whose Schwarzschild ra- 
dius is TO, and the scale factor a(t) is a solution of the 
Friedmann equation without the mass. It can be imme- 
diately seen that for a vanishing mass (H)) is reduced to 
the FRW metric, whereas when there is no cosmological 
source, i.e., a{t) = 1, we recover the Schwarzschild metric 
in isotropic coordinates. Describing an FRW spacetime 
with a massive source located at fixed comoving coordi- 
nates Xq, ([T]) fits well for the analysis of gravitational 
lensing where we, an observer at the coordinate origin, 
watch the deflection of light by a mass, i.e., a lens, mov- 



ing away from us according to Hubble's law. Then de- 
pending on the integrability of the geodesic equations we 
can perform step 2 either analytically or numerically. 

As a clearcut demonstration of our point, we will 
consider an exactly solvable case of a cosmological 
constant (A) driven universe, so that the scale factor is 
a{t) = e^* with H = -y/A/S. As one can guess nat- 
urally, a spacetime with a mass and A may be able to 
be described by the Schwarzschild-de Sitter(SdS) metric. 
Indeed, we can transform ([1]) into SdS by first transform- 
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and then redefining t hy t = t + f(f) with 

d/ Hr 

dr (1 _ m/f - H^f^)^l -m/f 

to reach a familiar form 

ds2 = -(l---H^f^]dP 
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But in the real world, the lcns(L), the source(S) and 
the observer(O) are all moving according to Hubble's law, 
and with ([5]) it may not be clear how to impose such a re- 
quirement, causing interpretational confusion. Therefore 
to facilitate the intuitive description of the system we go 
back to ID) and define "physical" spatial coordinates by 
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because the actual distance is (scale factor) x (comoving 
distance). Also, in weak gravitational lensing situations, 
TO, is much smaller than any other length scale under 
consideration, and we will work only up to 0{m). Then, 
the background of our interest is 
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FIG. 1: Lensing schematics. 
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where we align the coordinates to put the lens on the a;-axis. Note that the cosmic expansion is implemented in the 
coordinates themselves. Especially, the Hubble motion of L relative to O can be seen directly from the metric: It is 
moving along x = (— e^*g, 0, 0), where g is a constant related to the location of the lens at a certain moment of time. 
Gravitational lensing in SdS spacetime is not a new subject, and it was shown that A does not show up explicitly in 
observable quantities [3]: For a thin lens, the bending of the photon trajectory occurs for a period much shorter than 
H~^ ~ A"-'^/^, so that the spatial expansion due to A cannot have any effect practically. It does affect the propagation 
of the photon between S and L, and L and O, but the angular diameter distance properly absorbs it. But recently 
[3|-[J] took a fresh look at this old subject, claiming the appearance of 0(A) correction, and there have been prosQ 
and cons[6[ to their claim. Now we will revisit this problem with the metric of ([7]). 



II. LENSING IN A COSMOLOGICAL CONSTANT DRIVEN UNIVERSE 



We consider the motion of a photon in the spacetime of ([7]) (FIG. [T]). Corresponding Christoffel symbols can be 
calculated straightforwardly. Because of the symmetry of the system, we can confine the photon on the x-y plane and 
set z = 0. Also, the null condition 



0- -{l-^) + {l+'^){i^' - Hxf + iy' - Hyf}. 
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with a prime denoting a i-derivative and R — \J [x + ge-^*)^ + y^, can be used to simplify geodesic equations. To 
reduce the number of equations to solve, we parametrize the geodesies by t. Since t is not an affine parameter, our 
geodesic equations have the following form: 



x^" + VipX^^'xP' = kxf" . 



With the help of ([5]), the t-geodesic equation determines k to be 



k = H + —{-H{x^ + y2) + xx' + yy' + e"'q{x' ~ Hx)} , 
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and then x- and y-equations become 

= a;" - Hx + ^iHx'{x^ + y^) _ (^xx' + yy')x' + {xy' - x'y)y' + e"^q{Hx ~ x'){He"^q + 2x') + e"^q\ , 
= y" - Hy' + ^{Hy'ix^ + y') - (•^•^' + yy')y' - i^y' - x'y)x' + e"'q{Hy - y'){He"'q + 2x')] . (11) 

Next, we split the photon trajectory into the zeroth order piece(without L) and the first order one(with L, of (^(rn)): 
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x = xoA xi , y ^yo-i yi , (12) 
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where r is the typical length scale of the lens system and r ^ m. The zeroth order solution is trivial: 

xo = ai + aae^* , yo ^ h + feje^* . (13) 

In the absence of the lens, the boundary conditions are such that at i = ig the photon was at the location of S, 
{~xs,ys)i and at t — to it arrives at O at the origin, (0,0). Then, ai, 02, 61, 62 are fixed to be 

e"*°xs xs e"*oys ys . . 

Plugging p^ and (fT^ into pT|) . we get the first order geodesic equations: 
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whose solutions are 
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+ \/«? + ^?\/(«i + e^*(«2 + q)? + (&1 + e"'b2Y )) , 
yi = d2 + rfie^* - ''^°' + '^^ ^(fli + eg*(a2 + q)f + (61 + e^*&2)^ (16) 
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+ \/«? + 6?v^(ai + eff*(a2 + g))^ + (61 + 6^*62)^ )) 



When L is present, in order for a photon emitted at the source location to hit O at the origin, it should follow a 
bent trajectory, and the travel time is longer than io ~ ^s by an amount of 0{m). Let us say it reaches the origin 
at i = to + —ti- Then, we need to fix six unknowns, Ci, C2, rfi, ^21 9 and ti, to completely determine the photon 
trajectory. 

To be a trajectory of a photon (fT2| must satisfy the null condition: Plugging the solutions (fT3|) and ([16]) into ([8]) 
gives 

777- 

= H\al + bl)-l + -{Ai(i)(ai62 - a2&i) + A2{t){H\al + bj) - 1) + A3{t){aiC2 + 61^2)} , (17) 

r 

where A^'s are some complicated functions of t. Then, at 0{m'^), we get 

H^{al + bl) = l ^ e-^(*°-*^' = 1 - iJrs , (18) 



which relates the source-observer time difference, to — ts, to the source-observer distance, rs — v^^l + yf, in the 
absence of the lens. Next at 0{m) we have 

aiC2 + bid2 = 0, (19) 

and this is our first constraint on the unknowns. Note that (|18p can be generalized by considering ([5]) at 0{'m'^): 

d-X 
l + Hx=— =» H{to-t,)^logil + Hxo)-\og{l + Hx,) = -log{l + Hx,), (20) 

where x, is the location of the photon at t = t, . The boundary conditions in the presence of the lens give 4 constraints: 

Tji Tn vn 

x{ts) = -Xs = -Xs + fxi{ts) , x{to H ti) = = xo(io)— ti + —xi{to) , 
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The last piece of information necessary for connecting our algebra to the observation is the location of L: it was 
located at — tl at time t — tj^. In order for us to detect it at i = ioi '"l ^md ^l must satisfy (j20p . i.e., 
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With (HH), HH]) and ^^, solving ^^ and ^Bl gives the rest of the unknowns: 
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With ((24)) and (|25p . we can see that xi, yi and ii have an overall factor of r, so that the arbitrary parameter r does 
not appear in the full solution, ti has a special implication: — ii is related to the time delay due to the lens. 

Now that the photon trajectory is determined, we can calculate the angular location of the image of S observed by 



us: 
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where 13 — tan ^ -^ is the undeflected image location, and xl = 7'l/(1 ^ Hri^). This is the exact result up to 0{m). 
On the other hand the conventional Icnsing analysis (FIG. [2]), e.g., § IV of [7J, gives 
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FIG. 2: Lensing schematics in conventional analysis 



where d.'s are angular diameter distances defined by d, — a{t,)p, with p, being the corresponding comoving distance. 
Obviously 
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To compare (P7|) to our result (P5|) , we first note that our x^ , j/g and rg are distances at t = is and tl at t = f l . Then 
from ([5]), we immediately see that these are the same as the angular diameter distances, i.e., rs = dg, ri^ — di,, and 
therefore f3c = (3, 9c = 0- Now for small H, ([26]) becomes 



= /3 + -T7^|a;s - rfL + i?dL(a;s - ^l) + H^dl(xs - di^) + 0{H^) + 0{I3^)] + Oim") 



Using 
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we can rewrite (l29l) as 



where 



= a;s - rfL + ifdL(xs - dL) + i^'^Ll^s - ^l) + 0{H^) + 0{l3^) 



= P+^^^{l + O{H') + O{f3^))+O{m'). 
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Therefore, our result is in contradiction to the recent claims by ,3;l| which assert that there should be a 0{A) ^ 0{H^) 
correction to the conventional lensing analysis. 



III. DISCUSSION 



Through rigorous, close-to-the-reality but laborious 
derivations, we confirm that the conventional gravita- 
tional lensing analysis, although it may look loosely con- 
structed, is actually accurate, at least for thin, weak 
lenses in a A driven universe. [6|, which incorporated 
the "reality" in their analysis by employing FRW coor- 



dinates, has the same conclusion as ours. 

The differences between our work and 0-[j](and Q) 
can originate from the following: In 3|, 

1. the setup is not describing the actual observation: 
In the SdS universe, L and O must be moving 
relative to each other due to Hubble expansion, 
whereas they got an angle measured by an observer 
with fixed coordinates. 



2. their result is not written in terms of angular di- 
ameter distances, which are necessary to compare 
with the conventional results. 

In their follow-up paper i 4] , some effort was made to take 
these into consideration. But 

• their method to resolve problem 2, i.e., employ- 
ing the Einstein-Strauss scheme or the lens equa- 
tion, does not seem to be used properly: They are 
necessary in the Schwarzschild lensing because the 
Hubble expansion, which is not considered in the 
calculation of the bending angle, has to be imple- 
mented. On the other hand, the SdS spacetime is 
a complete playground for gravitational lensing be- 
cause the metric is equivalent to the McVittie solu- 
tion which already knows about both the lens and 
FRW. Then appropriate specifications of O, L and 
S should be enough for the full description of phe- 
nomena. Obviously, using Einstein-Strauss scheme 
or the lens equation together with SdS background 
is redundant. 

• the problem 1 was not addressed: Their angle ■0(in 
the first paper) or a(in the second paper) is mea- 
sured by a static observer(O). The instantaneous 



speed u of a geodesically moving observer(0/) can 
be calculated to be 



rL 



H^rl, 



(33) 



where tl is the radial coordinate of 0(0/) at the 
moment of the observation. Then, by relativistic 
aberration, the angle observed by 0/ is 
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assuming '0(V'') is small. To get the final result, 
we should convert quantities appearing in p4|) into 
angular diameter distances. This last step could be 
tricky, which makes the intuitive setup presented 
in the previous section favorable. 

The same rigorous procedure can be applied to generic 
FRW systems. For example, in the case of matter dom- 
ination with a{t) ^ t"^'^, once finishing tremendously 
more complicated intermediate steps, we can write the 
final result as 
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After relating dsh to other variables in a similar way as 
dSni), we get 

e = P+^^{l + 0{H'o) + Oif3^))+0{m^), (36) 

with Ho — 2/(3io), which still ascertains the accuracy 
of the conventional analysis. 

In addition to confirming the conventional results, our 
work can provide the exact formulae for errors coming 
from dropping higher order terms of 0{Hq) and O{0^), 
which might have importance in testing theories of long- 
range modification of gravity: Since our results are ob- 



tained under the assumption that the general relativ- 
ity(GR) is valid up to any scale, if precise measurements 
of high redshift lensing systems differ from our predic- 
tion, it may indicate that GR needs some modifications 
at large distances. At present, these theoretical errors 
are well buried beneath the measurement ones, but as 
our observational tools develop they may become signif- 
icant in the future. 
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